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Abstract 

In this review, we propose a summary of the most recent advances in the first-principles study of 
ferroelectric oxide epitaxial thin films and multilayers. We discuss in detail the key roles of mechan- 
ical and electrical boundary conditions, providing to the reader the basic background for a simple 
and intuitive understanding of the evolution of the ferroelectric properties in many nanostructures. 
Going further we also highlight promising new avenues and future challenges within this exciting 
field or researches. 
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I. INTRODUCTION 



During the recent years, numerous studies have been reported in the field of ferroelec- 
tric oxide nanostructures. These were motivated by the perspective of new technological 
applications and made possible thanks to combined spectacular developments of experimen- 
tal and theoretical techniques.- On the experimental side, it is now possible to grow oxide 
nanostructures with a control at the atomic scale and to measure their properties using local 
probes. On the theoretical side, different significant advances have been realized in the field 
of first-principles calculations that, combined to the recurrent increase of computational 
power, allow not only to reproduce accurately experimental measurements but also to make 
trustable predictions that are verified a posteriori. 

The study of ferroelectric nanostructures is nowadays a broad field of researches and it 
would be impossible to provide an exhaustive description of earlier works in a short review 
paper. Our purpose here is to provide a simple and comprehensive description of the origin 
of ferroelectric finite size effects and to highlight some recent results, focusing essentially on 
epitaxial thin films and superlattices. For a more complete and detailed description of recent 
first-principles achievements concerning ferroelectrics and their nanostructures, we refer the 
reader to recent book chapters^'^'^ and topical review papers.-'^>^'^'^>^ 

In the present review, we focus our attention on the infiuence of mechanical and elec- 
trical boundary conditions. Although other effects might also be important (chemistry of 
interfaces, defects ...), both effects play a key role and already allow to understand the main 
ferroelectric finite size effects when properly taken into account. The roles of mechanical 
and electrical boundary conditions are introduced separately and their respective infiuence 
is first illustrated in the paradigmatic case of a BaTiOs (the standard ferroelectric perovskite 
oxide material, see Sec. [Tll) thin film under realistic conditions, including its epitaxial growth 
on SrTiOs (Sec. IIII Al) . and the screening of its polarization charge by realistic electrodes 
(Sec. IIII Bl ). Once the basic physical effects have been introduced, the discussion is enlarged 
to other perovskite oxides epitaxial films (Sec. IIII CI) and superlattices (Sec. HVT ). Finally, 
recent works on other types of ferroelectric oxide nanostructures, such as nanoparticles or 
nanowires, are briefiy reviewed (Sec. [V|l . 
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II. BACKGROUND 



Barium titanate (BaTiOs) is a prototypical ferroelectric oxide.— At high temperature, it 
is paraelectric and crystallizes in the high-symmetry cubic perovskite structure as illustrated 
in Fig. [Ih: Ba atom is located at the corner of the cubic unit-cell while the Ti atom is at the 
center and is surrounded by an octahedra of O atoms, themselves located at the center of 
each face of the cube. On cooling, BaTiOs undergoes a sequence of three ferroelectric phase 
transitions to structures successively of tetragonal (T^ ~ 130°C), orthorhombic (Tc ~ 5°C) 
and rhombohedral (Tc ~ — 90°C) symmetry. In the following discussion, we will focus on 
the tetragonal phase that is stable at room temperature. 

The phase transition from cubic to tetragonal is characterized by the opposite shift of Ti 
and O atoms with respect to Ba, taken as reference. This polar atomic distortion is accom- 
panied with a small relaxation of the unit cell and yields a stable spontaneous polarization 
of 26 /iC/cm^ at room temperature. In the tetragonal phase, the cubic symmetry is broken, 
resulting not in a unique polar structure but in six symmetrically equivalent variants with 
polarization along the [100], [010] and [001] directions. The two variants along [001] ("up" 
and "down" states) are represented in Fig. [Th. What makes BaTiOs ferroelectric is not only 
the existence of energetically equivalent polar states but also the fact that it is possible to 
switch from one to another of these states by applying an electric field larger than the co- 
ercive field. Related to this switchability, the relationship between polarization and electric 
field is hysteretic in ferroelectrics, as illustrated in Fig. [T]:. 

In what follows, the pattern of cooperative atomic displacements associated to the phase 
transition will be referred to as ^, and the paraelectric, up and down states respectively 
will be labeled as ^ equal to 0, 1 and -1 respectively. For continuous evolution of ^, the 
internal energy U of the crystal has a typical double-well shape illustrated in Fig. [T]b. In 
BaTiOs, the well depth is typically of the order of 30 meV/cell. While hysteresis loops (Fig. 
d):) constitutes the measurement of choice to demonstrate experimentally ferroelectricity, 
the existence of a double-well shape for the energy is usually considered as the theoretical 
fingerprint of the ferroelectric instability. Making a Taylor expansion of U in terms of ^ 
around the paraelectric state taken as reference, and since the odd terms are forbidden by 
symmetry, we can write 
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FIG. 1: (a) Crystal structure of BaTiOs in its high temperature paraelectric cubic perovskite 
structure and in its room temperature tetragonal structure (for "up" and "down" polarization states), 
(b) Typical double- well shape for the internal energy of BaTiOs in terms of ^ ^ee text), fc) 
Hysteretic behavior of the polarization-electric field curve. Figure taken from Ref. 



= Ae. + Bit + ce +Dii+..., (1) 

where for the sake of simplicity we have assumed that atomic displacements in the tetragonal 
phase can take place only along z (condition that will be held all along this review unless 
otherwise stated), and we have cut the expansion in the eigth power of the atomic pattern of 
displacement. U presents a typical double- well shape with a negative curvature at the origin 
so that, in the previous equation, ferroelectricity can be associated to a negative value of A 
(A < 0). We note that, in the limit of a homogeneous polarization, the previous expansion 
in terms of microscopic atomic variable ^ is similar in spirit to that of Devonshire-Ginzburg- 
Landau (DGL) phenomenological theorj*^*^ in terms of the macroscopic polarization, since 

^ {i/n,)zu,, (2) 

where VLq is the unit cell volume and Z*^ is the diagonal component of Z* , the effective 
charge tensor associated to ^, so that part of the forecoming discussion can be transposed 
in that alternative context. 
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At the bulk level, a coherent microscopic model explaining the origin of ferroelectricity 
was first reported by Cochran,— who identified an intimate link between the structural phase 
transition and the lattice dynamics and introduced the concept of soft mode. Specifically, 
the underlying idea is that the lowest-frequency zone-center polar phonon mode in the 
paraelectric phase becomes softer as a function of decreasing temperature, and finally goes 
down to zero frequency, freezing in below Tc to generate the ferroelectric crystal structure. 
Cochran's theory of displacive phase transition was exhibited in the context of a shell-model 
and the softening of a given polar mode was explained from a competition between short- 
range and long-range Coulomb interactions. 

Most of these ideas were further confirmed from first-principles calculations^ performed 
within the density functional theory (DFT). An unstable polar mode (i.e. a mode with 
imaginary frequency u within the harmonic approximation) has been identified at F in the 
paraelectric phase of various ferroelectric oxides^ and, in prototypical cases like BaTiOa, a 
strong resemblance has been observed between the eigendisplacement vector of this unstable 
phonon and the ground-state distortion ^ (overlap of 99 % in BaTiOs or LiNbOa^). Since 
energy curvature is related to the square of the frequency, an unstable mode with imaginary 
frequency in the paraelectric phase is the indication of a negative curvature of the energy at 
the origin [A < in Eq. The amplitude of a; does not directly provide the double- well 
depth that additionally depends on anharmonic effects but measures, to some extent, the 
strength of the instability. The origin of the ferroelectric instability was also investigated 
at the first-principles level and the competition between short-range (SR) interactions and 
long-range (LR) Coulomb forces has been confirmed and quantified.— In this context, the 
anomalous values of the Born effective charges of ferroelectric oxides^*^ was shown as an 
essential feature to produce an unusually large destabilizing Coulomb interaction responsible 
for the ferroelectric instability. 

The balance of forces resulting in the ferroelectric instability is also well-known to be 
strongly sensitive to strain and external pressure. The sensitivity of the ferroelectric distor- 
tion to strain is inherent to the piezoelectric behavior and was highlighted from the early 
stage of first-principles computations.—"^ Samara et al?^ proposed that external hydro- 
static pressure reduces and eventually totally suppresses ferroelectricity since SR repulsions 
increases more rapidly than LR destabilizing forces under pressure. This latter result was 
recently reinvestigated at the first-principles level,— confirming a rapid disappearance of 
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ferroelectricity at moderate pressure and additionally pointing out an unexpected recovery 
of ferroelectricity at extremely high pressure. 

Such microscopic understanding of ferroelectricity at the bulk level allows to anticipate 
the existence of strong finite size effects in ferroelectrics since both SR and LR forces will 
be modified in nanostructures. The SR forces will be modified at surface and interfaces 
due to the change of atomic environment. The same is true for the Born effective charges 
that cannot keep their bulk value^^ as it was illustrated on BaTiOs thin films by Fu and 
coworkers.— Beyond that, the LR Coulomb interaction will be cut because of the finite size 
of the sample and will be strongly dependent on the electrical boundary conditions. Finally, 
since the ferroelectric instability is strongly sensitive to strains, it will be influenced by 
mechanical boundary conditions such as epitaxial strains.— 

These different factors can act independently to either enhance or suppress ferroelectricity. 
They will compete with each other in such a way that it is difficult to predict the ferroelectric 
properties of nanostructures without making explicit investigations. 

III. EPITAXIAL THIN FILMS 

During the recent years, most of the efforts devoted to understanding ferroelectric finite 
size effects concerned ferroelectric thin films. This was partly motivated by the observation 
at the end of the nineties of a polar ground-state in ultrathin films^i^^ with thicknesses 
well below 10-20 nm, that was still considered at that time as a realistic estimate of critical 
thickness for ferroelectricity.— 

It is now understood that the survival or disappearance of ferroelectricity in ultrathin 
films is not a truly intrinsic property: it strongly depends on the structure and chemistry 
of the interface between the ferroelectric and the substrate and/or the electrodes as well as 
on electrical (e.g. screening of the depolarizing field) and mechanical (e.g. epitaxial strains) 
boundary conditions. On top of that, additional features such as the finite conductivity of 
the film or the presence of structural defects and vacancies might also affect the ferroelectric 
properties. Although these latter features are potentially important, the proper treatment 
of their effect remains beyond the scope of first-principles simulations and is not yet fully 
understood so that these effects will not be further discussed here. For an overview of atomic 
relaxations and discussions of the stability of surfaces and interfaces we refer the reader to 
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previous reviews.-'^'^ In what follows, we will restrict the discussion on the role of mechanical 
and electrical the boundary conditions which, in most cases, provide a sufficient background 
to achieve relevant understanding of observed ferroelectric finite size effects. 

A. Mechanical boundary conditions : the epitaxial strain 

As mentioned in Sec. [Tll ferroelectricity in bulk compounds is well-known to be very 
sensitive to strains.—"^ In thin films, it is therefore expected that epitaxial strains will 
play an important role to monitor the ferroelectric properties. Coherently with that, strain 
engineering of the ferroelectric properties has been recently demonstrated experimentally: a 
polarization 250% higher than in bulk single crystals has been reported in BaTiOs epitaxial 
films,— and room temperature ferroelectricity has even been induced by strain in epitaxial 
SrTiOs thin films.^ 

In order to understand these effects, let us generalize our expansion of the internal energy 
in terms of additional strain Cij (where i and j are cartesian directions) degrees of freedom. 
In our paradigmatic example of a BaTiOa film epitaxially grown on a (001) cubic SrTiOa 
substrate we have mixed strain/stress boundary conditions: on the one hand the in-plane 
strains e^x = ^yy are fixed by the lattice mismatch between BaTiOs and SrTiOs, while 
Cxy = 0; on the other hand, the out-of-plane strain and the shear strains Cxz and Cyz are 
free to relax (condition of zero stress: azz = o'xz = CTyz = 0). Again, assuming for simplicity 
only an homogeneous polarization along 2;-direction, vanishing shear strains, and restricting 
the expansion to leading orders in ^ and e, the free energy functional to be minimized now 
reads 



[/(e, e) = Aez + Bit + Cil + Dil 

+ \Cii{'2el, + eL) + \Ci2{2el, + 4e,.e,,) 

+2gQexxil + i.9o + gi)(^zzil- (3) 

The terms in the first line correspond to the double- well energy of Eq. The terms in the 
second-line are the elastic energy while the terms in the third line arise from the coupling 
between ionic and strain degrees of freedom. They correspond to the so-called "polarization- 
strain coupling" and are at the origin of the piezoelectric response. It is clear from Eq. 
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([3]) that the polarization-strain coupling terms are responsible for a renormalization of the 
quadratic part of U that now takes the form 



[A + 2goe^^ + {go + Qi) e^J 



(4) 



Depending of the value of the parameters go and gi, and of e^x and e^z (deduced from 
the relation dU/dczz = which follows from the boundary condition a^z = 0), we see that, 
playing properly with the epitaxial strain conditions, it is possible to modulate the coefficient 
of in order to increase (resp. decrease) its negative value and therefore enhance (resp. 
suppress) the ferroelectric character of the film. 
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FIG. 2: Phase diagram of BaTiOs in terms of temperature and epitaxial misfit strain. The notation 
established by Pertsev et ai— is followed: p = paraelectric phase; r = monoclinic phase, related 
with the rhombohedral bulk phase; c = tetragonal phase with an out-of-plane polarization along the 
normal; aa = orthorrombic phase with the polarization along the cube face diagonal perpendicular 



to the normal. Figure taken from Ref. 
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The simplified free energy U in Eq. ([3]) can be generalized to include other orientations of 
^ as well as shear strains and the different parameters can be directly determined from DFT 
calculations at the bulk level.— The coefficients of the expansion for a series of ferroelectric 
perovskites have been reported by Dieguez et alM These energy parametrizations at K are 
the starting point for the development of model hamiltonians,-^^»^ that combined with clas- 
sical Monte Carlo simulations allow the first-principles based studies of temperature versus 
strain phase diagrams in perovskite oxide films, as illustrated for BaTiOa on a cubic (001) 
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substrate in Fig. [2l Although effects such as free top surface relaxation or interaction with 
the substrate are neglected, the phase sequence and the topology of the phase diagram are 
usually well reproduced within this approach. The predicted phase transitions temperatures 
are however typically too low compared to those experimentally measured. 

Alternatively, and in fact first applied to perovkskites by Pertsev and coworkers^ before 
first-principles investigations, similar phase diagram can be obtained from phenomenological 
Devonshire-Ginzburg-Landau (DGL) type expansion of the free energy in terms of polar- 
ization and strain. Results have for instance been reported for BaTiOa,— "^"^ PbTiOa,— 
SrTiOa,— "^"^ and Pb(Zr2;Tii_^)03 (PZT) solid solution.— This provides more direct access 
to the temperature behavior and gives very accurate results around the temperature/strain 
regime in which the model parameters were fitted (usually near the bulk ferroelectric tran- 
sition). In distant strain-temperatures regimes their predictions might be less accurate, 
however. For large strains or low temperatures the uncertainty increases and, as reported 
for BaTiOs, different sets of DGL parameters can provide phase diagrams qualitatively 
different far from the bulk transition region.—*^ 

As a general rule for the usual perovskites on a (001) substrate,— sufficiently large epi- 
taxial compressive (resp. tensile) strains will favor a ferroelectric c-phase (resp. aa-phase) 
with out-of-plane (resp. in-plane) polarization while distinct behaviors are predicted for the 
different compounds in the intermediate regime. Of course, there might be some discrep- 
ancies with this rule, such as the one pointed by Catalan and coworkers^ on PbTiOs films 
under tensile strains. 

For the BaTiOa film on a (001) SrTiOa substrate previously considered, the large com- 
pressive epitaxial strain (-2.38 %) will strongly increase the ferroelectric phase transition 
temperature and favor a tetragonal c-phase with an out-of-plane polarization. As further 
estimated by Neaton and Rabe,— such epitaxial contraint should enhance the polarization 
at K by more than 50 %. The sensitivity of polarization to strain is however not expected 
to be always so large: in highly polar Pb-based systems such as Pb(Zr2:Tii_2:)03 and PbTiOa 
where the ferroelectric displacements are already large, it was reported to be much more 
reduced than in BaTiOs, ~ similarly to what was also discussed for BiFeOs and LiNbOs.— 



In Ref. 



26l . like in many studies, the phase diagrams were obtained restricting the in- 



vestigation to the effect of the epitaxial strain on the ferroelectric degree of freedom ^, or 
equivalently P. In systems where different kinds of instabilites coexist (such as ferroelectric 
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and antiferrodistortive instabilities in SrTiOs,— or the interplay between in-plane ferroelec- 
tricity, antiferroelectric and antiferrodistortive distortions in PbTiOs free standing slabs^) 
the epitaxial strain can modify the competition between these instabilities so that all of them 
must be explicitely considered and generate much more complicated phase diagrams.-® Also, 
up to now, DFT calculations were restricted to mono-domain configurations while much com- 
plex domain structures can appear as highlighted using a DGL approacb^"^ or an effective 
Hamiltonian approach^"^ as it will be further illustrated in Section IIII Ci 

The effects imposed by the substrate on the lattice constant and geometry of the thin films 
can be maintained only up to a critical thickness, above which the elastic deformation 
can not be maintained any longer. Above he the strain energy is relaxed by the formation of 
misfit dislocations at the film-substrate interface. An estimate of he can be obtained using 
the Matthews-Blakeslee formula.— A priori, this structural limit imposes some restriction 
on the potential enhancement of the ferroelectric properties by the mechanical boundary 
conditions. However, it has been found that a coherent strain state can be preserved well 
above the theoretical critical thicknesses expected for strain relaxation.— Moreover, in some 
cases, it was reported that the ferroelectric transition temperature still has a value much 
larger than in bulk, even after the relaxation of the strain state.— 

In all previous approaches, only the presence of an homogeneous strain was considered. 
However, "flexoelectric" coupling of a strain gradient to the polarization,-^^^ and the related 
coupling of strain to a polarization gradient is receiving increasing attention. Although 
the characteristic scale of the flexoelectric coefficient is extremely small, it might affect the 
ferroelectric transitions in thin fllms. The subtleties of deflning flexoelectric coefficients as a 
bulk property^'^'^'^ and its computation in a first-principles framework is a challenge for 
theoreticians. 

A more exhaustive discussion of strain effects can be found in a recent topical review by 
K. M. Rabe.-^ 

B. Electrical boundary conditions : the depolarizing field 

As it was explained in Sec. [Tll ferroelectricity in perovskite oxides at the bulk level is 
a collective effect resulting from a delicate balance between SR interactions, that favour a 
paraelectric phase, and LR dipole-dipole electrostatic interactions, that favour a ferroelectric 
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distortion. So, it is expected that cutting these LR interactions by reducing the number of 
dipoles to interact with in a finite system, will aff'ect the ferroelectric behavior of the material. 

If we try to induce a homogeneous out-of-plane polarization in a (001) free standing 
BaTiOs slab in vacuum, preservation of the normal component of the electric displacement 
field at the interface will induce a depolarizing field £d = —P/eo, where eo is the permit- 
tivity of vacuum (unless otherwise stated we will use throughout this review the SI system 
of units). This depolarizing field totally suppresses the possibility of stabilizing an homoge- 
neous monodomain configuration with an out-of-plane polarization. 

Many applications rely on films with a uniform switchable out-of-plane polarization. So, 
after the previous discussion, it is clear that, to stabilize such a configuration, we need 
to screen the depolarizing field Sd- This can be achieved through two basic mechanisms: 
the first one is the compensation by "free" charges provided either by a metallic electrodes, 
semiconducting substrates, or even by ionic adsorbates^"^ (including in this category the 
screening by external applied fields that might be generated by the free charges or dipole 
layers)?^ the second one, and the only one if the substrate is an insulator,— is the breaking 
up of the system into 180° domains. Most of the fully first-principles calculations have been 
devoted to the screening by free charges, especially due to the very demanding computational 
effort that an explicit treatment of the domains implies, so that only the first issue will be 
discussed in more details here. 

In an early first-principles paper, Ghosez and Rabe,— using a model Hamiltonian ap- 
proach, showed that ferroelectricity might be preserved in ferroelectric slabs if Sd is per- 
fectly compensated. This was further confirmed by Meyer and Vanderbilt, at least for some 
terminations of the free standing slabs, in a fully first-principles atomistic calculation.— 

In principle, this total compensation of the depolarizing field could be achieved introduc- 
ing the ferroelectric film between two metallic electrodes in short-circuit. However, following 
an idea introduced in the early seventies by Batra,— "^"^ MethaA^ and coworkers in the con- 
text of a phenomenological model (see Fig. [3]), Junquera and Ghosez^ demonstrated from 
first-principles that, even in that case, compensation of the depolarizing field is incomplete, 
due to the finite screening-length at the metal/ferroelectric interface. 

Within this model, if a BaTiOs film gets a polarization normal to the ferroelec- 
tric/electrode interface, then the discontinuity of the polarization at the interface gives rise 
to a net interface charge density apor^ (Fig- EH) ■ In order to compensate for this polarization 
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FIG. 3: A schematic diagram of (a) a short-circuited electrode-ferroelectric structure with the 
spontaneous polarization displayed; (b) the screening charge distribution in the presence of perfect 
electrodes; its (c) charge distribution, (d) voltage and (e) field profiles in the presence of realistic 
electrodes. Please, note that here the film is taken to be a perfect insulator. Figure taken from Ref. 
71. 



charge density, a screening charge is induced in the electrode that spreads in practice over a 
finite distance from the interface (Fig. [3b), giving rise to an interface dipole moment density 
and a jump in the electrostatic potential at each ferroelectric/electrode interface (Fig. [3ld). 
This potential drop AV is shown to be proportional to the polarizatiop^ 



AV 



X 



(5) 



where Ag// is a constant of proportionality hereafter referred to as the effective screening 
length. Assuming a capacitor with two similar metal/ferroelectric interfaces, the total po- 
tential drop into the electrode/ferroelectric/electrode heterostructure amounts to 2AV. If 
the electrodes are in short-circuit (i.e. at the same potential), a residual field within the 
ferroelectric must compensate for this drop of potential (Fig. [3b), 



Sd 



2AV 

maQ 



(6) 



where mao is the film thickness in terms of the number m of unit cells of the ferroelectric 
perovskite oxide. Inserting now Eq. ([5]) into Eq. ^ 



2XeffP_ 

Eomao 



(7) 
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This shows that the field produced by the screening charge within the metal does not fully 
compensate for the bare depolarizing field £d = —P/sq: a residual depolarizing field remains 
that is proportional to P and to the ratio between the effective screening length and the 
film thickness. It will tend to zero for thick films (m oo) but can become sizable when 
the film thickness becomes comparable to Ag//. The ratio between the residual field into the 
ferroelectric (Eq. [7]) and the bare depolarizing field quantifies the quality of the screening 
and can be expressed as 1 — /3 = 2Ae///mao (as it will be further discussed in Sec. III-C). 

The depolarizing field can be computed averaging^'^'^ the microscopic electrostatic 
potential provided by a first-principles calculation for the whole capacitor. Details on the 
methodology and practical recipes to perform this nano-smoothing can be found in Ref. 



76l . Since the atomic positions, and therefore the induced polarization [Eq. ([2])], and the 
thickness of the film are input variables of the simulation, Xeff can be estimated using Eq. 
([7]) once £d is known. 

Junquera and Ghosez reported a value of Ae// = 0.23 A (one order of magnitude smaller 
than the interplanar distances in the heterostructure) for their frozen phonon calculations 
on a SrRuOs/BaTiOs/SrRuOs capacitor. It might be even slightly smaller when the atomic 
relaxation of the metal is properly treated.— It is worth noticing that the effective screening 
length introduced in this model is NOT an intrinsic property of the metal (like the Thomas- 
Fermi screening length), but a quantity that strongly depends also on the details of the 
ferroelectric/electrode interface such as the different chemical bonds formed at the junction, 
or the penetration of electronic wave functions of the metal into the ferroelectric and vice 
versa. As shown by Sai et al.J^ Xeff can vary and become very small for some specific 
metal/ferroelectric interfaces. Understanding and controlling this compensation mechanism 
of the polarization charges in order to stabilize a uniform switchable polarization in the film 
is one of the challenges in the field of thin film ferroelectrics.- 

Other authors22i2£ alternatively explain the origin of a residual depolarizing field inside 
the ferroelectric thin film from the presence of a "dead layer" at the metal/ferroelectric 
interface. Although evoking a different mechanism, the theory of the dead layer yields very 
similar mathematical expressions than the theory of the imperfect metal screening discussed 
above, the role of the effective screening length being played by the thickness of the dead 



layer (see for instance Ref. 
Levanyuk with those of Ref. 



2 
82 



or compare the formulas of Ref. 



81 



by Bratkovsky and 
by Chensky and Tarasenko). Stengel and Spaldin^ recently 
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reconciled these two apparently antagonist theories, providing a unified interpretation. Using 
a first-principles approach, they showed that the origin of an intrinsic dead layer with poor 
dielectric properties at the interface in oxide nanocapacitors is related to the hardening of 
the collective zone center polar modes that itself is produced by the incomplete screening of 
the depolarizing field. 

In the presence of a non-vanishing electric field, the proper energy to be considered 
and minimized is no more the internal energy U{^) of the crystal but a field-dependent 
energy potential F{^,£) that, on top of U{^), also includes an electrostatic energy term 
arising from the coupling between the polarization and the field.— This coupling term is 
proportional to —Sd ■ P- Since the depolarizing field is proportional to (— -P) [Eq. ([7])] the 
additional electrostatic term in F is proportional to and positive in sign. Remembering 
Eq. (I2l), we see that proper treatment of the electrostatic energy will produce an additional 
term in Eq. ([4]) of the form So, it clearly appears that, in a way similar to the 

epitaxial strain in Sec. IIII Al the incomplete screening of the depolarizing field is responsible 
for a renormalization of the quadratic coefficient of the energy. In this case however, the 
correction term is always positive (a > 0) and therefore its effect contributes to suppress 
ferroelectricity. 

The depolarizing field is inversely proportional to the thickness of the ferroelectric layer, 
Eq. (JT]). The thinnest the film, the larger the depolarizing field and the larger the energy 
penalty to stabilize a monodomain pattern of the polarization. In order to minimize the 
energy while remaining in a homogeneous configuration, the ferroelectric suffers a progressive 
reduction of the polarization that show itself well above the critical thickness. This point was 
first confirmed experimentally by Lichtensteiger et al^ for PbTiOs monodomain epitaxial 
films through X-ray diffraction (XRD) measurement of the evolution of tetragonality c/a 
with thickness. The tetragonality indeed provides indirect but easily accessible information 
on the polarization through the polarization-strain coupling [from Eq. (3), {c/a)p = (c/a)o + 
jP^] so that its measurement has become a standard technique to investigate size effects in 
thin films and superlattices.— "^"^ The reduction of the polarization was further confirmed 
by polarization data of Kim et a/.— in ultrathin SrRuOs/BaTiOs/SrRuOs capacitors free 
from passive layers,— and more recent X-ray photoelectron diffraction (XPD) experiments 
of Despont et al— 

As it was mentioned before, an alternative way of reducing the electrostatic energy is the 
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formation of 180° domains. In other words, going below the critical thickness predicted for 
monodomain configuration in ferroelectrics does not necessarily imply a paraelectric state.- 
The existence of 180° domains prevent the appearance of a net interfacial charge density.— 
This situation is well known in ferroelectric films on insulating substrates, where alternating 
180° stripe domains structures have been observed and characterized analyzing the satellite 
peaks in synchrotron X-ray diffraction measurements.—'^ Even on conducting substrates, 
where the free charges of the electrode should provide a substantial amount of screening, 
the formation of 180° domains has been proposed to explain concomitant suppression of 
polarization and recovery of full tetragonality in Pb(Zro.2Tio.8)03 on SrRuOs.— This is 
supported with recent works based on the Landau theory,— that conclude that domains 
always form in thin films, almost irrespective of the nature of the electrode and whether or 
not the screening carriers may be present in the ferroelectric itself. 

The situation is more complex experimentally, and whether the system breaks up into 
domains or not seems to be a very subtle issue. Lichtensteiger et al., using the same 
experimental setup, have observed how high-quality ultrathin films of PbTiOa grown on 
Nb-SrTiOs electrodes remain in a monodomain configuratiort^ (although with reduced po- 
larization and tetragonality) whereas they form domains when the electrode is replaced by 
Lao.erSro.ssMnOs/SrTiOs.— Even in the case of insulating substrates, the formation of do- 
mains depends on how quickly the sample is cooled. Monodomain phases of PbTiOs on 
insulating SrTiOs have been grown and the structure have been analyzed with a coher- 
ent Bragg rod analysis (COBRA) method,— when the samples were cooled from deposition 
temperature to room temperature over a period of 24 hours, in contrast to the 180° stripe 
periods observed before.—"^ Nevertheless, even in the latter case, the system can no more 
be properly called ferroelectric since, although locally polarized, its polarization is no more 
expected to be switchable.— 

C. First-principles modeling of ferroelectric capacitors considering mechanical 
and electrical boundary conditions. 

From the previous discussion, it is clear that appropriate treatments of the strain and 
of the incomplete screening of Sd are crucial to describe size evolution of the ferroelectric 
properties. Very recently, many DFT calculations as well as first-principles-based effective 
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Hamiltonian or atomistic shell-models simulations have been reported addressing the prob- 
lem of the thickness dependence of the ferroelectric properties in thin films, and the eventual 
existence of a critical thickness for ferroelectricity. 

A summary of various results obtained with DFT based methods is reported in Table-Hl 
In almost all these simulations, the prototype ferroelectric capacitor considered is made of 
an insulating layer of a ferroelectric perovskite oxide material (typically BaTiOs, PbTiOa 
or KNbOs) sandwiched between metallic electrodes (typically a conductive oxide such as 
SrRuOa or a noble metal such as Pt), although lately a free top surface and some adsorbed 
molecules have been also simulated.—"^ Typically the capacitor is built with two symmetric 
metal-ferroelectric interfaces, although very recently Gerra and coworkers^ have explored 
asymmetric heterostructures and linked the first-principles computations with phenomeno- 
logical approaches. The calculations were performed under short-circuit electrical boundary 
conditions for the electrodes (zero bias), a requirement that can be satisfied both in a fer- 
roelectric/metal superlattice geometry or in a periodic metal/ferroelectric/metal/vacuum 
geometryi^'i^ (it is only recently that a method has been proposed to address finite bias 
with periodic boundary conditions^). In all the previous approaches, a monodomain config- 
uration was considered, so the electrode was the only source of screening of the polarization 
charge. The mechanical boundary-conditions were imposed by fixing the in-plane lattice 
constant and geometry to the one imposed by the substrate. The final constraint depends 
on the author, but the most common approach was to assume a strain imposed by SrTiOa, 
the typical substrate on top of which these capacitors are epitaxially grown. The minimum 
of the Kohn-Sham energies were then searched either by a frozen-phonon method within the 
ferroelectric soft mode subspace (line minimization displacing the atoms in the ferroelectric 
layer by hand with an amplitude corresponding to a given percentage of the bulk soft mode 
as defined in Sec. [Tll) or by a full atomic relaxation of the atomic positions within the whole 
system. The results are summarized in Table-H] and deserve two important comments. 

First, although an exhaustive comparison of the results would certainly require to con- 
sider also other differences in the calculations, it is observed that DFT calculations performed 
within the Generalized Gradient Approximation (GGA) systematically overestimate the fer- 
roelectric character compared to those within the local density approximation (LDA): while 
LDA predicts the existence of a critical thickness for ferroelectricity that ranges between 
two and six unit cells (depending on the material, electrode, interface, and strain condi- 
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TABLE I: Review of the most recent fully first-principles simulations on the ferroelectric properties 
of ferroelectric ultrathin films in a monodomain configuration. Both computations within the local 
density approximation (LDA) and the generalized gradient approximation (GGA) to the density 
functional theory have been reported in the literature. CA functional refers to Ceperley-Aldei^ in 
the parametrization of Perdew-Zunger,— while PW91 and PBE stand for the GGA functionals of 
Perdew and Wan g^^'^^ and Perdew, Burke and Ernzerhof^^'^ respectively. Regarding the method, 

le S iesta code,— PW is the acronym 



1051, or DACAPOi^e in Refs 



60 



IE), 



NAO stands for numerical atomic orbitals as implemented in t 
for plane waves (as implemented in VASP-^*^ in Refs. [77 
and MBPP stands for mixed-basis pseudopotentials.-^ ay represents the in-plane lattice constant. 
If a theoretical in-plane lattice constant is used to fix the length of the in-plane lattice vectors of the 
heterostructure, it is computed within the same method and approximations than in the simulation 
of the whole heterostructure. Different methodologies yield to slightly different values of the lattice 
constant, tc stands for the critical thickness (if any) of the perovskite ultrathin film, in units of the 
number of cells of the ferroelectric perovskite oxide. 



Reference 


Heterostructure 


Method 


Functional 


Interface 


ay 


tc 


Junquera et al. [701 


SrRuOs /BaTiOs /SrRuOs 


NAO 


LDA (CA) 


SrO-Ti02 


3.874 


A (<^SrTi03) 


6 


Junquera et al. 


SrRuOs /PbTiOs /SrRuOs 


NAO 


LDA (CA) 


SrO-Ti02 


3.874 


A (f'-SrTiOs) 


6 


Gerra et al. [77] 


SrRuOs /BaTiOs /SrRuOs 


PW 


GGA (PW91) 


SrO-Ti02 


.3.94A(4h^;0^) 


3 


Umeno et al. [108] 


Pt/PbTiOs/Pt 


MBPP 


LDA (CA) 


Pt-PbO 
Pt-Ti02 


3.845 


A (o-SrTiOa) 


4 
6 








GGA (PW91) 


Pt-PbO 
Pt-Ti02 


3.905 


rt- V^pbTiOa-' 


No 
No 


Duan et al. [1051 


SrRuOa /KNbOa /SrRuOs 
Pt/KNbOa/Pt 


PW 


LDA (CA) 


SrO-Nb02 
Pt-Nb02 


3.905 


^ V"SrTi03^ 


4 
2 


Na Sai et al. [78.109[ 


SrRuOs/BaTiOs/SrRuOs 


PW 


GGA 


SrO-Ti02 


3.991 


A fo'^'P 1 
^ V"BaTi03^ 


> 4 










Ru02-Ba0 




> 4 




Pt/BaTiOa/Pt 






Pt-Ti02 
Pt-BaO 






> 4 

> 4 




SrRuOs /PbTiOs /SrRuOs 






SrO-Ti02 
Ru02-PbO 


3.905 


A 1 
rt- l"pbTi03^ 


No 
No 




Pt/PbTiOs/Pt 






Pt-Ti02 
Pt-PbO 






No 
No 


D. D. Fong et al. [60] 


SrRuOs/PbTiOs/vacuum 


PW 


GGA 


Sr0-Ti02 


"PbTiOs 


> 3 




SrRuOs/PbTiOs/OH, O or H 












< 3 




SrRuOs/PbTiOs/C02 












~ 3 




SrRuOs/PbTiOs/H20 












> 3 
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tions), calculations within the GGA often result in the absence of such a limit, especially 
in Pb-based perovskites, with ferroelectric ground states that are stable down to the ulti- 
mate thickness of one unit cell, eventually with an enhancement of the polarization.-^^»^ 
To interpret this trend it is worth to notice that, although GGA might a priori be consid- 
ered as an improvement over LDA, it does not necessarily produce more trustable results in 
the case of ferroelectric oxides. At the bulk level, the widely used GGA of Perdew, Burke 
and Ernzerhof (GGA- PBE)^ has been shown to strongly overestimate the ferroelectric 
character and, even, to yield an erroneous super-tetragonal structure in PbTiO^^'^'^ 
and BaTi03.~ Only the recently proposed improved GGA functional of Wu and Cohen^^ 
solves this problem and produces very accurate structures. The GGA results of Table-H] are 
not making use of this new functional so that their comparison with LDA results must be 
made with caution : indeed, the fact that GGA overestimates the ferroelectric character 
and concludes sometimes in the absence of critical thickness might be, at least partly, an 
artifact similar to that reported at the bulk level. 

Second, it is well-known that the bandgap of the Kohn-Sham particles as computed within 
DFT is significantly smaller than the experimental gap. For BaTiOs, typical LDA values 
are of the order of 1.6-1.9 eV— while the experimental value is of 3.2 eV. Although 
this DFT bandgap problem does not a priori constitute a failure of the theory that should 
internally correct for it, at least in "exact" DFT, in order to provide correct ground-state 
properties, it can yield unphysical results when working with approximate functionals. The 
band alignment at a typical metal/ferroelectric interface is illustrated in Fig. Hh where, for 
the purpose of the illustration, the Fermi level of the metal is assumed to be located roughly 
mid-gap, as typically expected from the experiment. The DFT bandgap underestimate 
prevents accurate predictions of the barriers simultaneously for electrons, and holes, 0p, 
and can in many cases produce pathological situations where the Fermi level of the metal 
is erroneously aligned with the conduction bands of the ferroelectric (0„<O), producing 
unphysical population of the conduction bands of the ferroelectric (Fig. Hb). Such a situation 



was, for instance, reported by the authors of Ref. 



83 for the Pt/SrTiO,-? interface, the 



structure of which had to be artificially modified to avoid the problem, or also ap pears 
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from the analysis of the partial density of states of the Pt/BaTiOa interface in Ref 
Moreover, even in "normal" cases where the DFT calculation reproduces a small barrier 
for electrons in a non-polarized configuration (Fig. Hb), the situation remains problematic 
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since the underestimate of the electron barrier artificially lowers the breakdown field: if 
during the relaxation of the ferroelectric state, the residual depolarizing field exceeds the 
erroneously low DFT breakdown field, a spurious transfer of electrons from the metal to the 
ferroelectric conduction bands will be observed at the interface, modifying its properties and 
producing an unphysical relaxed structure. For the SrRuOs/BaTiOs/SrRuOa nanocapacitor 



of Ref. 170|, we obtained within the LDA for the non-polar configuration (pp^"^ = 1. 4:7 eV and 
0^^^ = 0.11 eV, while the use of an empirical "scissors" correction (i.e. a rigid shift in energy 
of the BaTiOs conduction bands of 1.62 eV, to reproduce the experimental bandgap) provides 
a more realisitic value 0^"^ = 1.73 eV.— The artificially low value of 0^^^ prevented us 
to perform full atomic relaxations of the ferroelectric state but was not addressed by other 
authors reporting full relaxation on the same system. 

All these issues highlight that special care must be taken when dealing with 
metal/ferroelectric interfaces and invite to be specially critical when interpreting the re- 
sults. An important challenge at this point for the theoriticians is to identify an alternative 
functional avoiding the bandgap problem and its consequences on the band alignment, while 
preserving accurate description of the structural and ferroelectric properties. The Bl hybrid 
functional recently proposed by Bile et al.-^ might constitute an promising option. 

In spite of the technical difficulties highlighted above, first-principles DFT simulations 
remain the most accurate model currently available for the theoretical study of ferroelec- 
tric capacitors and provides a wealth of informations at the atomic level. Calculations are 
however performed at K, so that the study of the temperature dependence of the struc- 
tural, ferroelectric, and piezoelectric properties is not accessible. Also, DFT calculations 
are computationally very demanding and the effort required to deal with such complex het- 
erostructures, involving interfaces between dissimilar metal/insulator materials, is at the 
limit of present capabilities. To bypass these intrinsic limitations and address more com- 
plex questions, alternative less accurate methods have been proposed that remain based on 
first-principles results. These include effective Hamiltonian,-^^^ shell model,— i^*^ and 
molecular dynamics simulations combined with coarse-grained Monte Carlo simulations,— 
where the different parameters are directly fitted on DFT results. 

Wu and coworkers^"^ have performed effective Hamiltonian simulations on 
Pb(Zro.5Tio.5)03 free standing slabs. The choice of the Ti and Zr composition is interesting, 
since it is close to the morphotrophic phase boundary of bulk PZT. Since the calculations 
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FIG. 4: Schematic band alignment in a typical metal/ferroelectric/metal capacitor. In the para- 
electric state, due to the DFT bandgap problem, the Schottky barrier </)„ is underestimated (b), and 
can eventually produce pathological situations where (pn < (c) . In the ferroelectric state (P / 0), 
the DFT bandgap problem will additionally be responsible for a lowering of the breakdown field £c 
to S^^'^ ■ If the residual depolarizing field Sd associated to the ground-state polarization exceeds 
S^^"^ , a spurious transfer of electrons from the metal to the ferroelectric conduction bands will be 
observed at the interface and produce unphysical results. 

were done in the absence of surface charge screening, films thicker than 4 unit cell form 
180° domains in the direction normal to the surface to screen the depolarizing field, so the 
average component of the out-of-plane local mode vanishes, independently of the strain. 
The depolarizing field does not influence the in-plane polarization versus strain behaviour. 
A rich strain-temperature phase diagram, including a monoclinic and a stripe domain phase 
not present in the bulk parent material was predicted. 

For the same composition of Pb(Zro.5Tio.5)03, Kornev et al.^"^- extended the effective 
Hamiltonian in order to incorporate explicitly the presence of a free surface and the effects 



of a partial compensation of the depolarizing fleld. As in Ref. [92, in this first attempt. 



only the net component of the out-of-plane polarization at the surface layers was considered 
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to compute the homogeneous depolarizing field and the corresponding contributions to the 
energy. The treatment of Sd was then refined by Ponomareva and coworkers,— who cor- 
rected the dipole-dipole interactions of the bulk model hamiltonian by summing explicitely 
contributions of the local dipoles present within the film, and derived in this way an exact 
expression for the depolarizing energy and field at an atomistic level in any low dimensional 
ferroelectric structure under open circuit (OC) boundary conditions. If the polarization lies 
along a non-periodic direction, the OC boundary conditions naturally lead to the existence 
of a maximum of the depolarizing field that is later screened by hand. The amount of 
screening and therefore the amplitude of the residual field can be tuned with a parameter, 
P (see Section III-B), that plays a role similar to the effective screening length in the model 
described above and ranges between P = (OC boundary conditions) and /? = 1 (perfect 
short circuit boundary conditions). This model has been applied to different kinds of low 
dimensional ferroelectrics.— 

the interplay between the screening of the depolarizing field and the strain 



In Ref. 
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imposed by the substrate was also explored. In the limit of perfect screening, a uniform 
polarization along z is predicted, while the behaviour of the in-plane components of the 
polarization depends on the strain, following the general rule described in Sec. IIII A[ a 
compressive strain supresses the in-plane polarization while a tensile strain favours the pres- 
ence of a aa-phase. Therefore, the phase change from a monoclinic phase (the so-called 
Ma phase with {Px) = {Py) ^ 0, {Pz) > {Px)) for sufficiently large tensile strains to a 
tetragonal phase with (P^) = (Py) = 0, {Pz) > for stress free or compressive strains. 
On the other extreme case of open-circuit boundary conditions, the average polarization 
along z vanishes independently of strain, while the in-plane polarization aligns with the 
[010] direction {{Py) ^ 0, {P^) = 0) for the stress free, displays a monoclinic Mc phase 
{{Px) 7^ {Py) 7^ 0, {Py) > {Px)) under tensile strain, or vanishes for compressive strain. 
For intermediate values of the depolarizing field, the polarization continuously rotates and 
passes through low-symmetry phases. Interestingly, these phases are stable only for compo- 
sitions lying near the morphotrophic phase boundary of bulk PZT, disappearing in favour of 
a tetragonal phase for high enough Ti-concentrations. The fact that the average out-of-plane 
polarization vanishes under open circuit boundary conditions does not mean that the local 
dipoles along z vanishes. On the contrary, they might have a large value, close to the bulk 
one, but forming domains. The process of formation of these domains was monitored with 
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respect the screening of the depolarizing field. When 8d is large enough, some "bubbles" 
(cylindric nanodomains having local dipoles that are aligned in an opposite direction with 
respect the one in the monodomain configuration) nucleate and grow laterally till 180° stripe 
domains are formed. 

Model hamiltonian techniques have the additional advantage of a trivial treatment of 
external electric fields. Starting from a five unit cell thick thin film geometry under com- 
pressive strain, where c-oriented 180° stripe domains are stabilized even assuming a realistic 
amount of screening of the depolarizing field, the evolution of the domain structure under an 
electric field has been studied for Pb(Zro. 5X10.5)03^ and BaTiOs,— giving some light into 
the mechanisms to saturate the films. The path is essentially the same for both perovskites 
oxides and consist of three different steps: {i) the growth of the domain parallel to the field 
at the expense of the other up to a first threshold value of the external field; {ii) the for- 
mation of the nanobubbles defined before; and {iii) the contraction of the nanobubbles till 
the final collapse into a monodomain configuration at a second critical value of the external 
field. The main differences between BaTiOs and Pb(Zro.5Tio.5)03 are {i) the orientation of 
the domain walls: the stripes alternate along [110] for BaTiOs rather than along [100] for 
Pb(Zro.5Tio.5)03, and {ii) the "hardness" of the polarizarion in BaTiOs to rotate (BaTiOs 
films profoundly dislike significantly rotating and in-plane dipoles). 

The formation and response of the "bubbles" to external electric fields and the recent 
predictions of domains of closure^'^'^ (very common in ferromagnets but detected in 
ferroelectric thin films only recently^), have bridged a fascinating parallelism between fer- 
roelectric and ferromagnetic domains. Although the shape and size of the "bubbles" (ellipti- 
cal in ferroelectrics and spherical in ferromagnets), the abruptness of the domain walls (one 
lattice constant thick in ferroelectrics and much more gradual, over many atomic planes, in 
ferromagnets), or the exceptionally small periodicity and constancy with the field of the fer- 
roelectric domain periods differ from the counterpart ferromagnetic properties, they might 
not be so different as traditionally considered, despite the profound differences between elec- 
trostatic and magnetostatic interactions. Domains of closure have been predicted using a 
first-principles effective hamiltonian for Pb(Zro.4Tio.6)03^ asymmetrically screened (grown 
on a nonconducting substrate and with a metal with a dead layer as top electrode), and using 
a Landau-Ginzburg phenomenological approach for a PbTiOs thin film,— both asymmet- 
rically and symmetrically coated with insulating SrTiOs. These predictions have obtained 
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further credit after full first-principles simulations on SrRuOs/BaTiOs/SrRuOs ferroelectric 
capacitors by Aguado-Puente and Junquera,— where the domains of closure are obtained 
even for a symmetrical metal/ferroelectric/metal capacitor, where the metallic plates should 
provide significant screening. As it happens with ferromagnetic domains, ferroelectric do- 
mains follow the Kittel law (that states that the domain width is directly proportional to 
the square root of the sample's thickness), at least for thicknesses thicker than 1.2 nm.— 

The effective Hamiltonian approach is now considered as a standard method and has 
been used to address various other questions. The phase diagram of PZT films as a func- 
tion of temperat ure a nd Ti composition under open-circuit boundary conditions has been 
reported in Ref. Il38l . Other systems like BaTiOs have also been considered for which the 
strain/temperature phase diagram has been obtained and the effect of incomplete screen- 
ing of £d has been discussed.^ Although most studies focused on [001] epitaxial films, the 
infiuence of the growth direction on the properties was also investigated.— 

Besides all the previous theoretical frameworks, molecular dynamics combined with 
coarse-grained Monte Carlo simulations have been recently used to study a very important 
tecnological problem: the swithching of the ferroelectric polarization.— These simulations 
show that the shape of the critical nucleus is square and not a thin triangular plate as sug- 
gested by the prevailing Miller- Weinrich model. The critical nuclei grow in an anysotropic 
way, with activation barriers for sideways and forward growth much smaller than for nucle- 
ation. The values for the activation energies and activation fields predicted by the molecular 
dynamics simulations are one order of magnitude smaller than those expected from the 
Miller- Weinrich theory, mainly due to the new slanted interface model for the surface of 
the nuclei, and the domain velocities are in very good agreement with the ones measured 
experimentally.— 



IV. SUPERLATTICES 



Artificial ferroelectric superlattices are nowadays considered as an interesting alternative 
to ferroelectric thin films^ since, as it is discussed below, they allow fine tuning of the ferro- 
electric properties while maintaining prefect crystal ordering. The large variety of materials 
that can be combined (ferroelectrics, incipient ferroelectrics or regular insulators) and the 
way these materials can be ordered within the structure (bicolor and tricolor superlattices 
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or even more complex structures with composition gradients) also offers tremedous scope 
for creating artificial ferroelectric materials with possibly new and tailor-made properties. 

During the recent years, short-period epitaxial multilayers com- 
bining a ferroelectric and an incipient ferroelectric material (like 
BaTiOs/SrTiOs 4i-42 .i42.i43.i44.i45.i46. 147. 148.i49.i50. 151.152.153 PbTiOa/SrTiOa 86.87 .154.155 

KNbOs/KTaOs,— '^'^'^ and related tricolor systems like 

BaTiOa/SrTiOs/CaTiO^'i^'^) have retained most of the attention both at the 
theoretical and experimental level. In these systems, proper handling of the mechanical and 
electrical boundary conditions appears crucial to monitor the ferroelectric properties in a 
way very similar to what was previously discussed for thin films. 
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FIG. 5: Layer by layer polarization profile along a ferroelectric BaTiOs/SrTiOs/CaTiOs multilayer 
as obtained from DFT first-principles calcualtions. A polarization similar to that of BaTiOs is 



induced in the non-ferroelectric SrTiOs and CaTiOs layers. Figure taken from Ref. 
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First, the alternation of ultrathin layers is particularly suited to impose coherent epi- 
taxial strains and achieve "strain engineering" of the properties. The compressive epitax- 
ial strain produced by a SrTiOs substrate on BaTiOs layers contributes to explain po- 
larization enhancement with respect to the bulk BaTiOs value in BaTiOs/SrTiOs bicolor 
superlattices^'^'^ and even in tricolor— superlattices, despite the fact that part of the 
system is not nominally ferroelectric. The epitaxial strain additionally produces significant 
shift of the phase transition temperature as for instance reported for BaTiOs/SrTiO^i^ and 
PbTiOs/SrTiOs^ superlattices. Ferroelectricity is even expected in SrZrOs/SrTiOs super- 
lattices on SrTiOs, in spite of the fact that neither SrZrOs nor SrTiOs is ferroelectric at 
the bulk level.— Although it was not stated by the authors, this can probably be partly 
explained from the large compressive strain imposed on the SrZrOs layers by the SrTiOa 
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substrate. 

Second, electrostatic coupling between the layers is also playing a key role. On the 
one hand, polarizing the ferroelectric BaTiOs layers only in a BaTiOs/SrTiOs superlattice 
would produce a large polarization gradient at the interfaces, inducing huge depolarizing 
fields and electrostatic energy. On the other hand, polarizing SrTiOs also has an energy 
cost since this material is not spontaneously polarized. The ground-state of the system will 
therefore depend on the competition between different energy terms and can strongly evolves 
depending of the compounds involved and the superlattice period. 

In short-period superlattices where the ferroelectric alternates with a highly polarizable 
compound such as an incipient ferroelectric (like SrTiOs or KTaOs), a polarization will be 
induced in the latter in order to avoid depolarizing fields and the superlattice will typi- 
cally exhibit a homogeneously polarized configuration (see Fig. [5]). First pointed out by 
Neaton and Rabe on BaTiOs/SrTiOs superlattices, such a behavior was also theoretically 
predicted for PbTiOs/SrTiOs,— or even BaTiOs/SrTiOs/CaTiOgr^'^ short-period super- 
lattices. In many theoretical studies, the layer polarization was estimated roughly, multi- 
plying the atomic displacements into the layer by the bulk Born effective charges. A more 
rigorous and accurate approach has been proposed recently based on Wannier functions.— 

In a recent work, Dawber et al.— have demonstrated that key ferroelectric properties can 
be tuned over a very wide range in PbTiOs/SrTiOs superlattices: it was shown that the 
polarization can be adjusted from to 60 /iC/cm^ and the transition temperature from room 
temperature to 700°C by properly adjusting the volume fraction of PbTiOs. Moreover, a 
simple Landau-based model was proposed that allows to determine a priori the properties 
of the system, so opening the way to the production of samples with ferroelectric properties 
designed for particular applications. It was also reported that the presence of SrTiOs between 
the PbTiOs layers in the superlattice improve the electrical properties and allows to avoid 
the large leakage currents typically reported in PbTiOs thin films while maintaining perfect 
crystal structure. 

The simple arguments described above have been shown to provide reasonable descrip- 
tion of (PbTiOs)n/(SrTiOs)s superlattices over a relatively wide range of PbTiOs thick- 
nesses (for n ranging from 1 to 53).— They are however expected to apply mainly to rel- 
atively short-period systems. In longer period superlattices, other behaviors are expected. 
In KNbOs/KTaOsT^ it was theoretically shown using a shell-model that the electrostatic 
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coupling between layers (responsible for a homogeneous polarization throughout the system) 
will progressively decrease as the period increases, yielding a progressive reduction of KTaOa 
polarization. In BaTiOa/SrTiOa,— effective Hamiltonain simulations have confirmed that 
short-period superlattices behave like a homogeneous-like material with a phase diagram 
that resembles that of a (001) BaTiOs thin film under short-circuit conditions while inho- 
mogeneous atomic features and original domain patterns were predicted in longer period 
systems. Appearance of domain structures in ferroelectric/paraelectric superlattices was 
also discussed by Stephanovich et al.— 

The previous discussion implicitly focused on the evolution of the out-of-plane polar- 
ization. In some cases, an in-plane component of the polarization can also eventually ap- 
pear. This was predicted in KNbOs/KTaOs superlattices in which the epitaxial compressive 
strain is not strong enough to force KNbOs (with a rhombohedral ground-state) to become 
tetragonal.— In a way compatible with a symmetry lowering observed experimentally,— an 
in-plane polarization component was also predicted in the SrTiOs layers of BaTiOs/SrTiOs 
superlattices under in-plane expansion.— The main difference between the out-of-plane and 
in-plane polarization components is that while the out-of-plane components of the different 
layers are strongly electrostatically coupled as previously discussed, the in-plane components 
behaves much more independently and remain confined in the individual layers in which they 
are induced by the epitaxial strain conditions imposed by the substrate. 

The interest for ferroelectric superlattices is not limited to the fact that they allow the 
tuning of the ferroelectric properties as discussed above but also results from their poten- 
tiality to generate totally new behaviors or phenomena. In Pb(Sco.5+xNbo.5-a;)03 (PSN)^ 
and Pb(Zri_a;Ti2:)03 (PZT)^^ superlattices, it is expected that symmetry lowering can yield 
new stable polar phases and strong enhancement of the functional properties. Unusual ther- 
moelectric properties and nonergodicity were also predicted in PZT superlattices.— In very 
short-period PbTiOs/SrTiOs superlattices an unexpected recovery of ferroelectricity has 
been observed in the limit of ultrathin PbTiOs layers that cannot be explained from the 
simple arguments reported above.— "^i^. Beyond such bicolor systems, tricolor superlat- 
tices, which break the inversion symmetry^ are also expected to potetially exhibit unusual 
and attractive properties. 

First-principles based simulations contributed to make significant progresses in the un- 
derstanding of the behavior of ferroelectric superlattices and even opened the door to a 
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possible theoretical design of artificial ferroelectric nanostructures with predetermined and 
tailor-made properties.— 

V. NANOPARTICLES AND NANOWIRES 



Although most of the recent experimental and theoretical efforts on ferroelectric nanos- 
tructures concerned thin films geometries, the growth and characterization of ferroelectric 
nanowires,— nanotubes,— nanocolumns,— nanoislands,— or even "nanosolenoids"— 
are at an incipient and exciting state. The geometries of these nanostructures favor novel 
inhomogeneous polarization and strain configurations, strongly dependent of the electrical 
and mechanical boundary conditions. 

DFT first-principles studies of such low dimensional structures are, however, compu- 
tationally very demanding. Other recent studies made use of first-principles based shell- 
models^ and model hamiltonian methods.— For instance, domain structures, polarization 
and coercive fields of nanoscopic particles of BaTiOa have been carried out using such 
methodologies. 
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FIG. 6: "Vortex-like" structure of the local polarization in a (12 x 12 x 12) BaTiOs quantum dot. 
The figure corresponds to a cut along the y = 6*'* plane. Although the macroscopic polarization 
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letoroid moment of the vortex structure has a well-defined finite value. Figure taken 
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Perhaps, the most remarkable result is the prediction of novel exotic geometric ordering of 
ferroelectricity in BaTiOgr^ and Pb(Zro.5,Tio.5)0^^'^ nanoscale disks and rods, where the 
the local polarization rotate from cell to cell, forming a "vortex-like" pattern. Although the 
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macroscopic polarization of these quantum dots vanishes at zero apphed field, they exhibit 
a new order parameter, the toroid moment of polarization, introduced by Naumov et alr^ 
and defined as 

i 

where pi is the local dipole of cell i located at and N is the number of unit cells in the 
simulation. Since the local dipoles in the vortex can rotate either clockwise or anticlockwise, 
the toroid moment might adopt two different values pointing in opposite directions, so that 
one bit of information might be stored by assigning one value of the Boolean algebra ("1" 
or "0") to each of these states. The direction of the toroidal moment can be efficiently 
controlled by a transverse inhomogeneous electric field generated by two opposite charges 
located away from the studied dot,— opening exciting oportunities for nanomemory devices. 
Technologically, the ferroelectric vortex phase promises to increase the storage density of 
ferroelectric random access memories (Fe-RAM) by five orders of magnitude. 

In order to investigate the possibility of engineering the dipole patterns and vortex struc- 
tures, a study of the phase diagram as a function of temperature of these ferroelectic nan- 
odots when they are embedded in a polarizable medium was reported by Prosandeev and 
Bellaiche.— Six different phases were found, depending on the ferroelectric strength of the 
material constituting the dot and of the system forming the medium. Some of the novel 
phases are remarkable, for instance the coexistence of a toroidal moment and a spontaneous 
polarization at low T when a soft ferroelectric dot is immersed in a medium that is ferroelet- 
rically harder than the dot. Medium driven interactions between dots were demonstrated 
with adjacent vortices rotating in an opposite fashion ("antiferrotoroidic" phase). 

The phase transition from the vortex structure (with the local dipoles circularly ordered) 
to a ferroelectric phase (with all the local dipoles pointing along the same direction and 
with an out-of-plane polarization) in a Pb(Zro.5Tio.5)03 nanocylinder by an electric field 
perpendicular to the plane of the vortex has been recently addressed.— The transition 
between the two phases is through an intermediate state that generates a lateral toroid 
moment and breaks the macroscopic cylindrical symmetry. The observation of this new phase 
yielded to the discovery of a stricking collective phenomenon, that leads to the annihilation 
of the ferroelectric vortex in a peculiar azimuthal mode. 
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A first step to explain the existence of this toroidal moment in ferroelectric nanoparticles 
from a DGL theory has been taken by Wang and Zhang^ A more exhaustive review of all 
these computations can be found in Ref. Q. 

In between the two-dimensional (2D) thin films and superlattices discussed in Sees. [1111 
and HVl and the zero-dimensional (OD) nanoparticles described above, the case of one- 
dimensional (ID) nanowires was also recently addressed. Model hamiltonian simulations 
on Pb(Zr,Ti)03r^ and full first-principles simulations on BaTiOs nanowires^ have demon- 
strated how long-range ordering exists in these ID systems even at finite temperature and 
under zero external field. These results challenged traditional statistical models that claimed 
that ID lattices with particle-particle interactions decaying faster than r~^,n > 3, are im- 
possible to have phase transitions at finite T. As it happened in thin films and nanoparti- 
cles, the mechanical and electrical boundary conditions are essential in the nanowires. The 
presence of uncoated free surfaces yields to the existence of depolarizing fields along two 
directions and forces the collective polarization to lye along the longitudinal direction. A 
critical diamenter for the polarization was estimated around 1.2 nm for BaTiO^^ and 2.0 
nm for Pb(Zr,Ti)03,— below which the ferroelectricity was suppresed. At smaller radii, low 
atomic coordinations at the surface produce a contraction of the unit cell that is responsible 
for the suppression of the ferroelectric distortion at the critical radius. Nevertheless, even 
below this threshold size of the diameter, ferroelectricity might be recovered playing with 
the mechanical boundary conditions, by applying appropriate tensile strain. 

Previous first-principles simulations on nanowires also contributed to get a new insight 
into the concept of the ferroelectric correlation volume, defined as the smallest volume within 
which atomic displacements must be correlated in order to decrease the energy of the crystal 
and produce a stable polar entity. Analyzing the critical diameter of the wire as a function of 
its length, the computations prove that is the anisotropic shape of the polar domain rather 
than its volume that determines the stability of its ferroelectric state.— 

VI. CONCLUSIONS AND PERSPECTIVES 

As illustrated throughout this review, various progresses have been reported during the 
last years concerning the description of the ferroelectric properties in oxide ultrathin films 
and superlattices. In this context, first-principles theory and modeling constituted an ef- 
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ficient guide and provided a significant support to experimental works. On the one hand, 
proper understanding of the crucial role played by the electrical and mechanical allowed to 
realize significant advances as discussed above. On the other hand, however, many issues 
still remain unclear so that this field of research is still at an exciting and incipient stage. 

In ferroelectric thin films, despite recent achievements, some important questions remain 
without clear answer and interesting perspectives ask for further confirmations. Only re- 
cently addressed from first-principles,— the mechanism and the quality of the screening at 
free surfaces with air remains unclear. Also the reason why some ferroelectric films on 
metallic electrodes stay monodomain with a reduced polarization while others break into 
domains is particularly intriguing.— The roles of defects, such as oxygen vacancies, and 
of leakage currents on the amplitude of the depolarizing field was up to now neglected in 
theoretical studies while they might play a role, at least in some cases. The mechanism of 
the switching and domain-wall motion in ultrathin films is an important question that has 
been studied in detail experimentally^'^'^'^ but that has not been addressed yet at the 
first-principles level, although the first-steps on multi-scale computations based on molecular 
dynamics and coarse-grained Monte Carlo simulations have been recently taken.— Recent 
studies reported interesting perspectives for ferroelectric tunnel junctions^'^'^ but no 
first-principles quantitative estimate of the tunneling current and of the electroresistance 
has been reported yet. Finally and without being exhaustive, the case of a ferroelectric 
layer between magnetic electrodes was proposed as a mechanism to achieve ferroelectric 
control of magnetisnt^i^ but is still widely unexplored. 

Ferroelectric superlattices constitute a promising alternative avenue that offers practical 
solutions to bypass some of the drawbacks of epitaxial ultrathin films (huge depolarizing 
fields, high leakage currents, strain relaxation), while preserving perfect crystal ordering 
and allowing to monitor the ferroelectric properties. Up to now, fine tuning of the polar- 
ization and phase transition temperature has, for instance, been illustrated playing with 
the superlattice period and composition but the most exciting perspectives probably rely 
on the potentiality of such systems to generate totally new and unexpected phenomena. 
Also, superlattices might offer particularly attractive opportunities in the emergent field of 
magneto-electric multiferroics.— 

The perspectives are numerous and first-principles simulations will certainly continue to 
play an active role in the field of nanoscale ferroelectrics in the future. In some cases, this will 
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mandatory require new theoretical developments such as, for instance, the identification of 
alternative and more accurate functionals or of new efficient approaches to compute complex 
properties. In order cases, the computer power will constitute the main limitation. In all 
cases, pursuing the constructive dialog initiated during the recent years between theorists 
and experimentalists will certainly appear as a major motivation. 
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